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By first-principles calculation based on the density funetional theory (DFT) with the modified 
Becke-Johnson loeal density approximation plus Hubbard U (MBJLDA-l-t/), we studied the 
band struetures of the Yb ehaleogenides YbO and YbS under ambient and high pressures. It 
was revealed that both YbO and YbS have a trivial band topology under ambient pressure, and 
a nontrivial band topology under high pressure. The topologieal phase transition is redueed 
by the pressure-driven single-band inversion between 5d- and 4/-orbitals at the time-reversal 
invariant momentum (TRIM) point X. A bulk Dirae eone eoexisting with a pair of metallie 
surface states on the [001] surfaee determined by tight binding model ealeulation with a 
slab geometry also demonstrates the nontrivial band topology of YbO and YbS under high 
pressure. 


1. Introduction 

Because of the emergent novel surfaee state protected by nontrivial band topology, 
researeh on topologieal elassifieation for gapped systems is surging,In the family of 
symmetry-protected topologieal order, the topological insulator (TI) protected by time- 
reversal symmetry attraets a lots of researeh interest on fundamental theory and realization 
in concrete materials for potential industrial applieations.^^^ Not only conventional semieon- 
duetor materials with strong spin-orbital eoupling (SOC),^ but also Kondo insulator materials 
are attraeting inereasing researeh interest.^ In 2010, on the basis of the periodic Ander¬ 
son Kondo lattiee model, Dzero et al. proposed the topologieal elassifieation in a simple 
eubie Kondo insulator.^ Their work showed that, with the fluctuating valence of a rare-earth 
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element, the transition between topological trivial and nontrivial states can be realized in 
a simple cubic Kondo insulator prototyped using SmBg.*^ Both first-principles calculation 
and experiments suggest that the simple cubic SmBg with a mixed valence state is a can¬ 
didate topological Kondo insulator (TKI).‘^ PuBg is also predicted as candidate TKI by 
first-principles calculation. 

Very recently, Weng et al. have predicted that the simple cubic YbBg is a candidate topo¬ 
logical insulator, while YbBi 2 with a rock-salt structure is a potential candidate topological 
crystalline insulator with a nontrivial mirror Chem number protected by mirror symmetry.^^’ 
Angle-resolved photoemission spectroscopy (ARPES) also suggested YbBe as a candidate 
TKI.^^ In fact, there are hundreds of rare-earth compounds RX with a rock-salt structure, 
with R referring to the rare-earth Ce, Pr, Nd, Pm, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and 
Eu, and X standing for the pnictide atoms N, P, As, Sb, Bi and the chalcogenide atoms O, S, 
Se, Te, Po.^"^ There are many semiconductors or semimetals with different combinations of 
R and X. In particular, for the chalcogenides of Eu, Sm, and Yb,^^ a fluctuating valence state 
is accessible under ambient or high pressure.These rare-earth compounds providing a 
high probability of finding more TKIs, although there are few TKIs predicted in rare-earth 
pnictides and chalcogenides with a rock-salt structure at present. 

In this work, we employ first-principles calculation based on the density functional the¬ 
ory (DET) with the modified Becke-Johnson local density approximation plus Hubbard U 
(MBJEDA-l-t/) to study the band topology of the Yb monochalcogenides YbO and YbS with 
a rock-salt structure (shown in Eig. 1). Experimentally, both YbO and YbS are insulators 
with band gaps of 0.32 and 1.40 eV respectively, at ambient pressure,^^’^° and YbO is on 
the margin of a mixed valence state.Our MBJEDA-l-f/ calculations predict that, at ambient 
pressure, both YbO and YbS are trivial insulators with band gaps of about 0.24 and 1.26 
eV, respectively, and both of them will undergo a topological phase transition to a nontrivial 
topological metallic state under high pressure. The insulator-metal transition under high pres¬ 
sure is consistent with the experimental result.^^’^° The novelty of this work is that we also 
predict that this transition is a topological phase transition, i.e., the insulating and metallic 
phases have different band topologies. The mechanism for this topological phase transition is 
pressure-driven band inversion at the time-reversal invariant momentum (TRIM) point X be¬ 
tween the 4/-dominating band and the 5d conduction band owing to the increasing bandwidth 
of 5<i-derived states under high pressure. At ambient pressure, the Yb atom in both YbO and 
YbS is divalent because of the electronic correlation mimicked by Hubbard U, which local¬ 
izes the 4/-orbitals and suppresses the charge transfer between the 4/- and 5<i-orbitals. The 
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suppressed eharge transfer between the 4f- and 5(i-orbitals renders the band topology of YbO 
trivial. Under high pressure, the bandwidth of the JJ-orbital elearly inereases. Once the band 
minimum of the 5<i-orbital is lower than the band maximum of the 4/-orbital, charge transfer 
between the 4f- and 5<i-orbitals will take place. The pressure-induced charge transfer will 
lead to a mixed valence state in YbO and YbS, which is a necessary condition for 

2. Methods 

The DFT calculations employed the all-electron, full-potential linearized augmented 
plane wave (FPLAPW) method with MBJLDA implemented in the WIEN2K code.^^’^^ To 
account for the on-site 4/-electron correlations, we applied the MBJLDA -l- U scheme that 
incorporates the on-site Coulomb repulsion U and the Hund’s rule coupling strength Jh with 
an approximation correction for self-interaction correction (SIC) for the Yb 4/-orbitals.^^ The 
spin-orbit coupling was included in the self-consistent calculations, and the on-site Hund’s 
exchange parameter J// was set to zero. We used RMT'xKmax=9-0, muffin-tin radii of 2.50 a.u. 
for Yb, and 2.33 (2.05) a.u. for S (O), and a 30x30x30 k-point Monkhorst-Pack mesh. We 
used the experimental lattice parameters a=4.86 and 5.68 A for the lattice parameters of YbO 
and YbS under ambient pressure. Under high pressure, i.e. reduced lattice, we present the 
calculated band structures with a=4.70 and 5.10 A for YbO and YbS, respectively, where the 
(i-band will inverse with the top/-band. 

The basis set consisted of the Yb 4/, 5d, 6s, and 6p valence states, the Yb 5^ and 5p 
semicore states (treated in another energy window), and the O, S 3s semicore and 3p valence 
states. 

3. Results and Discussion 

Ambient pressure The calculated band structures of YbO and YbS at ambient pressure are 
shown in Fig. 2. Both YbO and YbS are predicted as good metals by MBJLDA calculation, 
and the band inversion between 5d- and 4/-orbitals exists in both band structures. However, 
YbO is a semiconductor with a 0.32 eV band gap, and YbS is an insulator with a 1.40 eV 
band gap as determined by experiment.^^’^° While the MBJLDA calculation is crucial for the 
noninteracting topological insulators,^® and describes accurately the correlation effects of the 
5d states, it can’t describe well those of the more localized 4/ states, whose treatment still 
requires an additional Hubbard U value. 

The MBJLDA-l-t/ calculations with U of ~ 3.0 eV render a band gap of ~0.24 eV for 
YbO, while YbS needs a larger Hubbard U =7.0 eV to obtain a band gap of 1.26 eV. Our band 




Fig. 1. (Color online) Crystal structure of YbO and its Brillouin zone (BZ). 


gap determined by MBJLDA+i7 ealeulation is close to the experimental gaps of 0.30 and 1.40 
eV for YbO and YbS, respectively.^^’A larger U is required to obtain the experimental gap 
for YbS, whose lattice parameter is larger than that of YbO. This variable U is consistent 
with previous work on Yb,^^ where a decreasing U is required under increasing pressure, i.e., 
reduced lattice parameter. The band structures of YbO and YbS determined by MBJLDA+t/ 
calculation are shown in Figs. 3(a) and 4(a), respectively, in which SOC split /-bands into 
two sub-bands consisted of/5/2 and/ 7/2 states, respectively. At ambient pressure, d- and/- 
bands are well separated and there is no hybridization between d- and/-orbitals. Because of 
the empty J-band, the Yb ion is divalent. The band topology is trivial since there is no band 
inversion in the entire Brillouin zone (BZ). 

High pressure Even though both YbO and YbS are regular band insulators at ambient 
pressure, we can still induce a topological phase transition in both of them by applying exter¬ 
nal pressure. Since there is no band inversion in both YbO and YbS because of the absence 
of d-f hybridization and valence fluctuation at ambient pressure, we can apply high pressure 
approach to extend the wave functions of 5d- and 4/-orbitals and realize the charge transfer 
between 5d- and 4/-orbitals. Once the band minimum of the conducting 5d band is below 
the band maximum of the/7/2 band, charge transfer between 5d- and 4/-orbitals will take 
place, which will cause the system to be in a mixed valence state. The calculated band struc¬ 
tures with the reduced lattice parameters a=4.70 A for YbO and 5.10 A for YbS are shown 
in Figs. 3(b) and 4(b), respectively. The calculated band structures show that both YbO and 



















Fig. 2. (Color online) Band structures of YbO with experimental lattice parameters a=4.86 A(a) and YbS 
with experimental lattice parameter a=5.68 A(b) obtained by MBJLDA calculation. The weights of 5d and 4/ 
electrons are represented by dark and gray, respectively. 


YbS become metallic without a global band gap under high pressure. In contrast to the band 
structures of YbO and YbS at ambient pressure, the bandwidth of 5(i-dominating bands in¬ 
creases significantly under high pressure, and band inversion indeed takes place around the 
X point. This band inversion will change the parity of occupied bands at the X point.'^*’ Since 
the 4/-orbital is odd under spatial inversion, while 5<i-orbital is even, the d-f hybridization 
vanishes at the X point, and each band at the X point has definite parity eigenvalues. The 
parity eigenvalues at all time-reversal invariant momentum points, including one T, three X 
points and four L points, are listed in Table I. Although both YbO and YbS are metals un¬ 
der high pressure, we still can define the Z 2 topological invariant owing to finite local band 
gap over the entire BZ.'*® In fact, except the well-known topological metal bismuth, several 
more topological metals have been predicted, such as the nonmagnetic GdBiPt,"^^ AuTlTe 2 , 
and SmS under high pressure.Since both YbO and YbS have crystal inversion symmetry, 
the Z 2 invariant can be calculated as the product of half of the parity (Kramers pairs have 
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Fig. 3. (Color online) Band structure of YbO with experimental lattice parameter a=4.86 (a) and 4.70 A(b) 
obtained by MBJLDA+t/ calculation, with U-3.Q eV. The weights of 5d and 4/ electrons are represented by 
dark and gray, respectively. 


identical parities) eigenvalues, Si, for all the occupied states at the TRIM points, 

8 

i=i 

v=0 reveals trivial band topology in insulating YbO and YbS at ambient pressure, while 
the nonzero v=l in the metallic phases of YbO and YbS reveals the nontrivial band topology 
of YbO (YbS) with lattice parameter a=4.70 A (5.10 A). With a shorter lattice parameter, the 
band inversion will involve a J-band and a lower/-band, but the number of band inversion at 
the X point still is one. The predicted metallic properties of YbO and YbS under high pres¬ 
sure are consistent with experimental results.We also note that there is a band crossing 
between the SJ-dominating band and the 4/-dominating band. However, the band crossing is 
absent in compressed SmS and SmO. This band crossing can be explained by the absence of 
long-range d-f hybridization. By the extended 3D Bernevig-Hughes-Zhang (BHZ) model on 
a face-centered cubic (FCC) lattice, the nearest-neighboring (NN) d-f hybridization allows 
two gapless nodes in the r(0,0,0)-Z(0,0,2;r) direction..Because of the shorter lattice param¬ 
eter of SmS, the next-nearest-neighboring (NNN) d-f hybridization cannot be ignored, and 
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Fig. 4. (Color online) Band structure of YbS with experimental lattice parameter a=5.68 (a) and 5.10 A(b) 
obtained by GGA+U calculation, with U-1.0 eV. The weights of 5d and 4/ electron are represented by dark 
and gray, respectively. 


Table I. The product of parity eigenvalues for all occupied states at the T, three X, and four L TRIM points 
in the BZ, and the total parity product of all TRIM points. The + and - denote even and odd parity, respectively. 


r 3X 4L Total 

insulating YbO -- ____ + 

insulating YbS -- ____ + 

metallic YbO - + + + ---- - 

metallic YbS - + + + ---- - 


this NNN hybridization will make the band structure locally gapped in entire BZ. In the YbO 
and YbS, the NNN d-f hybridization can be ignored because of the larger crystal lattice. 

Surface state By first-principles calculation using maximally localized Wannier functions 
(MLWFs),'^'*’"^^ we extract the on-site energy of/-orbitals and J-orbitals of Yb and the hopping 
parameters. By tight binding (TB) calculation, we exactly reproduce the band inversion at X 
point and the band crossing predicted by MBJLDA calculation. The [001] surface state is 
determined by TB calculation with a slab geometry consisting of 150 layers of YbS unit 
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Fig. 5. Surface state at (001) surface determined by TB calculation with slab geometry. 


cells. The surface state is shown in Fig. 5, and the 4/-eleetron poeket is slightly downward 
shifted to elearly show the surfaee state. We ean find a pair of metallie surfaee state around 
the r point. The bulk Dirae eone also is also present near the f point. This surfaee state is 
very similar to the [010] surfaee state of NasBi,'^^ in whieh a pair of metallie surfaee states 
eoexist with the bulk Dirae eone. 

4. Summary 

In summary, we prediet pressure-indueed topologieal phase transition in the ehaleo- 
genides YbO and YbS with a roek-salt strueture by DFT ealeulation, despite different mag¬ 
nitudes of eorrelation for 4/ eleetrons in YbO and YbO resulting from a lattiee differenee. At 
ambient pressure, beeause strong eorrelation loealizes the/-orbitals, eharge transfer between 
4/ and 5d orbitals is suppressed, and both YbO and YbS are regular band insulators. How¬ 
ever, under high pressure, the bandwidths of 4/ and 5d bands inerease, and 5d band is below 
one or several 4/ bands at the X point with inereasing pressure. The 4/ and 5d band inversion 
of YbO and YbS under pressure ehange the global band topology and render a topologieal 
metallie phase. Although the insulator-metal transition of YbO and YbS under high pressure 
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is well-known, we prediet this phase transition is also aoeompanied by a ehange in the global 
band topology. The bulk Dirae eone and a pair of metallic surface states coexist on the [001] 
surface state. 
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